Phase-diffusion dynamics in weakly coupled Bose-Einstein condensates 
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We study the phase-sensitivity of collisional phase-diffusion between weakly coupled Bose-Einstein 
condensates, using a semiclassical picture of the two-mode Bose-Hubbard model. When weak- 
coupling is allowed, zero relative phase locking is attained in the Josephson-Fock transition regime, 
whereas a tt relative phase is only locked in Rabi-Josephson point. Our analytic semiclassical 
estimates agree well with the numerical results. 



Bose-Einstein condensates (BECs) of dilute, weakly- 
interacting gases are currently used as a testing ground 
for the quantitative study of various condensed-matter 
models. In addition to offering a remarkable degree 
of controllability, they open the way for exploring non- 
equilibrium dynamics, far beyond small perturbations of 
the ground state. Whereas this regime is inaccessible 
in the equivalent condensed matter realizations, due to 
the high Fermi energy, highly excited states are naturally 
produced in BEC experiments and their dynamics can be 
traced with great precision. 

One such example is the 'phase-diffusion' [J, [3, |3, 0, H, 
0, 01 between two BECs prepared in a coherent state and 
consequently separated. Principle aspects of this process 
are captured by the two site Bose-Hubbard Hamiltonian 
(BHH). Defining di, a\ as bosonic annihilation and cre- 
ation operators for particles in condensate i = 1,2, the 
corresponding particle number operators are fii = a|di, 
and the BHH takes the form |l 



H = eL,- JL^ + ULi 



(1) 



where ~ {a\a2 + aldi)/2, Ly = {a\d2 — aldi)/{2i), 
and L2 = {fii — ri2)/2, generate an SU{2) Lie algebra. 
Bias potential, coupling, and collisional interaction ener- 
gies are denoted as e, J, and U, respectively. We have 
eliminated c-number terms that depend on the conserved 
total particle number N = hi-\-n2. Below we use for rep- 
resentation the Fock space basis states \£, m)^, which are 

the joint eigenstates of L and L^, with £ — N/2. 

Experimental procedures allow to prepare the sys- 
tem in any desired SU{2) coherent state \9,ip) — 
g-i(pi,^g-iei,„|^^£^^_ The 9 rotation is reahzed by strong 
coupling J S> NU , while the phase is tuned by switch- 
ing the bias e, thus inducing coherent phase oscillations 
^, 11]. In particular it is possible to prepare the equal- 
population coherent states 
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In phase-diffusion experiments, the preparation stage is 
followed by a sudden separation of the condensates, so 



as to obtain two equally populated modes of a symmet- 
ric (£=0) double well with a definite relative phase 
Due to the interaction terms, different Fock states 
oscillate with different frequencies, and the evolution un- 
der the Hamiltonian ([1]) with u = NU/J ^ 1 leads to 
the loss of single-particle coherence, quantified by the 
multi-shot mean fringe- visibility g[^2 — \{Lx)\/£- For 
fully separated condensates the fringe visibility decays 
as 5^2^^) = e-(*/*^)' with U = {UVI)-^ [H, i i fl. 
Long time revivals are obtained at tr = tt/U, as con- 
firmed experimentally This behavior is insensitive 
to the initial phase ip, and merely refiects the Binomial 
Gaussian-like distribution of the occupation, implied by 
Eq. 1^. 

In reality we still may have some inter-mode cou- 
pling J, which is small compared with the collisional en- 
ergy NU 0,13, Sim. Recent experiments in ID BECs, 
have explored the loss of single-particle coherence under 
such circumstances 0, resulting in the decay to a non- 
vanishing equilibrium value of g[^2 ■ note that since 
these experiments focus on phase fiuctuations associated 
with one dimensionality , their detailed description goes 
beyond the two-mode BHH 0, , and that dephasing 
can also be caused by thermal noise 0. Other exper- 
iments have found (p-sensitivity in the full merging of 
two separated condensates ll|, where a relative ip—ir 
phase led to significant heating losses compared to (p=0. 
The common expectation in these experiments is that 
small coupling between the condensates will lead to their 
phase-locking and suppress the loss of single-particle co- 
herence. 

Generally motivated by this experimental interest we 
here study the f dependence of the fringe visibility evo- 
lution when the two condensates remain weakly coupled 
during the hold time, assuming the prototype two-mode 
BHH modeling. In Fig. [1] we plot the numerically cal- 
culated g[^2 both the |7r/2,0) and the |7r/2,7r) prepa- 
rations. For the (f=0 preparation, the expected phase- 
locking is obtained even by a very weak coupling, leading 
to a nonvanishing asymptotic value of the fringe visibility, 
with some small oscillations. However, if the system is 
prepared in the (p=7r state, phase locking is not attained 
if u is moderately large, and the fringe visibility exhibits 
apparently complex dynamics with a somewhat different 
characteristic timescale. The corresponding short-time 
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FIG. 1: (Color online) Evolution of fringe visibility with 
A^=1000 particles, starting from the spin coherent states 
j7r/2,0) (a-d), and \tt/2,tt) (e-h). The coupling parameter 
is M = lO"* (a,e), 10^ (b,f) 10^ (c,g), and 10 (d,h). 
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FIG. 2: (Color online) Schematic drawing of the classical 
phase-space energy contours (a) and the corresponding spec- 
trum (b), for A'^ = 20 and u — 10. WKB energies (red x) 
are compared with exact eigenvalues (blue +). Dashed lines 
indicate slopes ujj for low energies, u>x for near-separatrix en- 
ergies, and UJ+ for high energies. Squares mark the eigenstates 
plotted in Fig. 3. 



dynamics was previously studied using truncated and lin- 
earized models [ij. Our purpose in the present work is 
to provide both qualitative and quantitative longer time 
analysis of this phase-sensitivity of the phase-diffusion 
process, using a semiclassical phase-space method. 

Interaction Regimes — It is instructive to use a 
Wigner function p( 6, (f) for the representation of the 
quantum states [12l |. The Wigner function that repre- 
sents a |0, Lp) coherent state resembles a minimal Gaus- 
sian centered at the corresponding angle on the BHH's 
spherical phase space. Using semiclassical reasoning one 
can argue that the Wigner functions of the BHH eigen- 
states (with e = 0) are concentrated along the contour 
lines of the Gross-Pitaevskii classical energy functional. 
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Furthermore, within the framework of the WKB approx- 
imation [isl . 14 1 the quantization of the energy is im- 
plied by the condition A{En) = (47r/A^)(n -I- 1/2) where 
A{E) is the phase-space area enclosed by a fixed energy E 
contour, and Att/N is the Planck cell. [Below we mea- 
sure A{E) in Planck cell units]. The phase space pic- 
ture implies that one has to distinguish between three 
regimes according to the value of the dimensionless in- 
teraction parameter u = NU / J [1, 0]: (a) The linear, 
weak-interaction Rabi regime u < 1; (b) the intermedi- 
ate Josephson regime 1 < u < N'^; and (c) the extreme 
strong-interaction focA; regime u > iV^. In the Josephson 
regime the spherical phase-space is split by a figure-eight 
separatrix trajectory (see Fig [2]), to a 'sea' of Rabi-like 
(blue) trajectories and two interaction-dominated nonlin- 
ear 'islands' (green). In the Fock regime the phase-space 
area of the sea becomes smaller than Planck cell, and 
therefore effectively disappears. 

Spectrum — The energy landscape for u 3> 1, as 



implied by Eq.Q, consists of sea levels that extend 
from the bottom energy — —£J, and of island lev- 
els that occur in almost degenerate pairs and extend 
up to the upper energy E+ « i'^U. The border be- 
tween the sea and the islands is the separatrix energy 
E^ ^iJ. The oscillation frequency uj{E) = [A'{E)]-^ 
around the minimum-energy (stable) fixed point (7r/2, 0) 
is the plasma / Josephson frequency 



ojj = oj{E^) = y/{J + NU)J « VNUJ (4) 

If we had u < I then the opposite phase-space fixed point 
(7r/2,7r) would be also stable with w+ = \/{J — NU)J. 
But for u > 1 it bifurcates, and replaced by an unsta- 
ble fixed point at the separatrix energy, accompanied by 
the two twin stable fixed points (arcsin(l/u), tt), located 
within the islands. For w 3> 1 these stationary points 
approach the poles (arcsin(l/M) = 0, tt) and the oscilla- 
tion frequency becomes ~ NU . The significance of 
the various frequencies is implied by the WKB quantiza- 
tion (Fig. 2). At low-energy we have a non-degenerate 
set of Josephson levels with spacing loj. Due to the non- 
linearity the spacing becomes smaller as one approaches 
the separatrix energy (see next paragraph), while the 
high-energy levels {Em = Um?) are doubly-degenerate 
with spacing 2U m that approaches w+ as m ^ ^. 

Nonlinearity — The motion in the vicinity of the 
separatrix is highly non-linear. This is implied by the 
non-linear variation of the phase space area near the sep- 
aratrix: 



\AiE) - A{E^ 



E-E^ 



log 



NJ 



E-E^ 



(5) 



We stress that Planck cell units are used for A{E). On 
the basis of this expression, the WKB quantization con- 
dition implies that the level spacing at the separatrix is 
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FIG. 3: (Color online) Wigner function of the separatrix state 
(a) and its n and tp projections for the eigenstates specified in 
Fig. 2 (b,c). Solid red lines depict the near separatrix eigen- 
state, while dash-dotted blue and dashed green lines corre- 
spond to the ground and most excited ('cat') states, respec- 
tively. Also shown in (d) are the relative-number variance and 
the fringe-visibility for ground state (blue dash-dotted line) 
and most excited state (dashed green line) with 7V=1000. 



finite and given by the expression 
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It is important to notice that in the strict classical limit 
this frequency becomes zero 15|, while for finite N it 



differs from the Josephson frequency only by a logarith- 
mic factor. The non-linearity can be characterized by a 
parameter a that reflects the E dependence of the os- 
cillation frequency uj{E), and hence upon WKB quan- 
tization it reflects the dependence of the level spacing 
En+i — E„ « u}{En) on the running index n. The im- 
plied definition of this parameter and its value in the 
vicinity of the separatrix are: 



a{E) = -'^^uo{EfA"{E) 
Lu an 
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(7) 



Note that this parameter is meaningful only deep in the 
Josephson regime, where it is less than unity. As shown 
below, the parameter a controls the amplitude of the 
fluctuations in the case of a ip=TT preparation. Strangely 
enough, in the classical limit {N oo with fixed u) 
the non-linear effect becomes more pronounced as far as 
Eq.([6l) is concerned, but weaker with regard to Eq.d?]). 

Eigenstates — (Fig. [3]) In the Rabi regime, the 
eigenstates of Hamiltonian approach the Lx eigenstates 
\£,m')x, with eigenvalues —Jm' proportional to the rela- 
tive number difference between the even and odd super- 
positions of the modes. In particular, the lowest eigen- 
state \i,i)x and the highest eigenstate \£,—i)x, are the 
coherent states \n/2,0) and |7r/2,7r), respectively, with 
binomial m distributions in the \£, m)z basis, approach- 
ing the normal (Gaussian/Poissonian) distribution for 
large i. As u increases beyond unity, a transition is 



made into the Josephson regime. The relative-number 
variance AL^ of the ground state decreases continuously 
in this regime from its coherent y/N /2 value. However, 
the coherence of t;he ground state remains close to 
unity and the relative-phase is well-defined throughout 
the Josephson regime, justifying the use of mean-field 
theory to depict the Josephson dynamics of ground state 
perturbations. The coherence of the ground state is only 
lost in the Fock regime, where the true ground state ap- 
proaches the relative site-number state \i,Q)z- 

(/3=0 preparation — The above discussion implies 
that in the case of a coherent preparation |7r/2,0), any 
weak coupling beyond J > U/N would suffice to lock the 
relative phase and to prevent the loss of single-particle co- 
herence [5| . Using a phase space language this means that 
the iy9=0 preparation is located at the minimum of the sea 
region, around a stable fixed point (7r/2,0) , and resem- 
bles the ground state of the Hamiltonian. The expected 
small-oscillation frequency for the evolution of this co- 
herent state is wj, and hence that of the visibility is 2ujj 
(see Fig. 4). 

(/3=7r preparation — For the coherent preparation 
|7r/2, tt) the picture is quite different. In the Rabi regime, 
this state coincides with the |€, —tjx most excited eigen- 
state, and hence it is possible to have phase locking 
around the stable fixed point (7r/2,7r), where the char- 
acteristic frequency is u;+. However, as soon as u > 1, 
this fixed point bifurcates, and the Wigner function of 
the eigenstates has to equally populate the two twin is- 
lands in phase space, hence forming a cat state with 
bimodal m distribution. This crossover is refiected (see 
Fig. [3|) by the sharp increase in the relative-number vari- 
ance to the cat-state value AL^ = A''/2, and by the 
equally sharp decrease of its single-particle coherence. 
The evolution of the coherent preparation |7r/2, tt) in the 
vicinity of the unstable point is further analyzed below. 

Dynamics — The phase-sensitivity of phase-diffusion 
originates from the instability of the hyperbolic fixed 
point (—1,0,0) of the classical phase space, as opposed 
to the stability of the elliptic fixed point (1, 0, 0) y. Here 
we would like to go beyond this heuristic observation 
and understand quantitatively the long time recurrences 
in Fig. [H The coherent state |7r/2,0) is a superposi- 
tion of a few low-energy states. Since the level-spacing 
is fixed on the Josephson frequency wj, the fringe vis- 
ibility carries out clean harmonic oscillations. By con- 
trast, the state |7r/2,7r) consists of levels around the sep- 
aratrix energy, which are not equally spaced, resulting 
in the quasi-periodic oscillations of Fig. [He)-(h). The 
observed timescale of these oscillations (Fig. 4) agrees 
well with Wx- For the purpose of quantitative analysis 
we have to expand the initial state if) in the En ba- 
sis. A reasonable estimate for the envelope function 
P{En) = \ {En\i})? = trace(p(")p('>)) can be obtained us- 
ing a semiclassical approximation. The Wigner function 
of the nth eigenstate is approximated by a microcanon- 
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FIG. 4: (Color online) Probability distribution P(-E„) for the 
preparations |7r/2,0) (a) and |7r/2,7r) (b), with A'^=1000 par- 
ticles, and u = 10"^. Numerical expansions (o) are compared 
with the semiclassical estimates of Eq.([8} and Eq.(|9]) (solid 
lines). In (c), the mean frequency and the long time aver- 
age of the (I/a:) oscillations in Fig. l(b,f) (A for |7r/2,0), □ 
for |7r/2, tt)) are compared to the semiclassical analysis of the 
dynamics in phase space (solid blue and dotted red lines, re- 
spectively). Error bars reflect the associated dispersion. 



ical distribution p(")((/7,6l) = uj{En)&(E(B,<^) - E„) and 
the coherent state is approximated by a minimal Gaus- 
sian. With appropriate approximations for u 3> 1 the 
phase space integration gives for the Lp—Q preparation 



P(-B) = 2Io 



2u 



-(2+l/2ti)_E 



(8) 



with E = (E — E^)/ J, and for the ip—T: preparation 

1 



TT V LVJ 



2+— I \E\ 

2u 



(9) 



with E = {E — Ey.)/J. In the above, Iq and Kq are 
the modified Bessel functions. As shown in Fig. [H^a) 
and Fig. [IJb) this analytical estimate agrees remarkably 
well with the exact numerical diagonalization, with the 
exception of the separatrix state whose phase-space dis- 
tribution (Fig. [3lja)) somewhat deviates from the micro- 
canonical ansatz. Due to the even parity of both prepa- 
rations (|7r/2,0) is always even and |7r/2,7r) has (—1)^ 
parity) , the occupation of odd n states vanishes and the 
occupation of the even n states is twice the semiclassical 
estimate. Using these expressions we can estimate the 
participation ratio of the preparation, which comes out 
M - NU/ujj « 

Fluctuations — Regarding ?/' as a superposition of 
M energy states we can estimate the fluctuations using 
the following reasoning: In the linear case, if the energy 
levels are equally spaced, there is only one M = 1 ba- 
sic frequency and the motion is strictly periodic. If the 
energies were quasi-random, then there are A4 = M—1 



independent frequencies, and the motion is quasiperiodic 
with relative variance ~ 1 /M . But what we have for the 
TT preparation are energy levels that are characterized by 
a nonlinearity parameter a. Then, if the nonlinearity 
parameter is not very small, the effective number of in- 
dependent frequencies is = ^/aM, which interpolates 
between the linear and the quasi-random estimates. As 
shown in Figs. 4(c) and 4(d), this analysis is consistent 
with our numerical simulation. 

Summary — The vast majority of recent BEG inter- 
ference experiments related to phase-diffusion, are car- 
ried out in the Josephson regime. Within it, the dynam- 
ics of single-particle coherence is expected to strongly 
depend on the initial relative-phase between the par- 
tially separated condensates. Zero relative phase leads 
to phase- locking even for u as large as N'^, whereas a tt 
relative phase preparation results in the loss of fringe- 
visibility all the way down to m ^ 1. Using a semi- 
classical phase-space picture we related this behavior to 
the classical phase-space structure and the ensuing WKB 
spectrum, obtaining analytical expressions for the enve- 
lope functions of these two coherent preparations in the 
eigenstate basis, for u in the Josephson regime. These 
expansions agree well with numerical calculations and 
enable the accurate prediction of the amplitude of fringe- 
visibility oscillations. 
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